Study of Loschmidt Echo for a qubit coupled to an XY-spin chain environment 
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We study the temporal evolution of a central spin-1/2 (qubit) coupled to the environment which is 
chosen to be a spin-1/2 transverse XY spin chain. We explore the entire phase diagram of the spin- 
Hamiltonian and investigate the behavior of Loschmidt echo(LE) close to critical and multicritical 
point(MCP). To achieve this, the qubit is coupled to the spin chain through the anisotropy term as 
well as one of the interaction terms. Our study reveals that the echo has a faster decay with the 
system size (in the short time limit ) close to a MCP and also the scaling obeyed by the quasiperiod 
of the collapse and revival of the LE is different in comparison to that close to a QCP. We also show 
that even when approached along the gapless critical line, the scaling of the LE is determined by 
the MCP where the energy gap shows a faster decay with the system size. This claim is verified by 
studying the short-time and also the collapse and revival behavior of the LE at a quasicritical point 
on the ferromagnetic side of the MCP. We also connect our observation to the decoherence of the 
central spin. 

PACS numbers: 05.50.+q,03.65.Ta,03.65.Yz,05.70.Jk 



I. INTRODUCTION 

Recent years have witnessed a tremendous progress in 
studies of quantum information theoretic measures [l|, l2[ 
close to a quantum critical point (QCP) [J-Q- Quantities 
like concurrence [fjHH , negativity [{§ LL1| , quantum fidelity 
pTl4l5| . quantum discord [16| etc., have been found to 
capture the ground state singularities associated with a 
quantum phase transition (QPT); for recent reviews see 

01. 

On the other hand, the studies of decoherence namely, 
the quantum-classical transition by a reduction from a 
pure state to a mixed state have also attracted the at- 
tention of physicists in recent years [19l - l22| . In this con- 
nection, the concept of Loschmidt echo (LE) has been 
proposed to describe the hypersensitivity of the time evo- 
lution of the system to the perturbation experienced by 
the environment to which it is coupled |23l - |27| . The mea- 
sure of the LE is the modulus of the overlap between two 
states that evolve from the same initial state \tf> ) under 
the influence of two Hamiltonians Hq and Hq + 5, where 
5 is a small perturbation, given by 



L{t) = \{i> \e i{Ho+5)t e- iHot \^)\ 2 - 

In some of the recent works, attempt has been made 
to connect these two fields by studying the behavior of 
the LE close to a QCP as a probe to detect the quantum 
criticality. Quan et al, studied the decay of LE using 
the central spin model where a central spin-1/2 (qubit) 
is coupled to the environment which is chosen to be a 
transverse Ising chain of N spins in such a way that it is 
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globally coupled to all the spins of the spin chain through 
the transverse field term [28[. The coupling to the qubit 
leads to the perturbation term S defined above and con- 
sequently, the time evolution of the spin chain initially 
prepared in its ground state, gets split in two branches 
both evolving with the transverse Ising Hamiltonian but 
with different value of the transverse field. This results 
in the decay in the LE. It has been observed that the LE 
shows a sharp decay in the vicinity of the quantum crit- 
ical point of the environmental spin chain; at the same 
time at the QCP, the LE shows collapse and revival as 
a function of time with the quasiperiod of revival of the 
LE being proportional to size of the surrounding. This 
study has been generalized to the case where the envi- 
ronment is chosen to be a transverse XY spin chain and 
the behavior of the LE has been studied close to the Ising 
critical point driven by the transverse field [29|, |3(| . 

Rossini et al |3l|, studied a generalized central spin 
model in which the qubit interacts with a single spin of 
the environmental transverse Ising spin chain and it has 
been shown that the decay of the LE at short time is 
given by the Gaussian form exp(— Ft 2 ) where the decay 
rate T depends on the symmetries of the phases around 
the critical point and the critical exponents. For instance, 
for such systems with local coupling, it has also been re- 
ported that r has a singularity in its first derivative as 
a function of the transverse field at the QCP [3l|. In a 
subsequent work [32(, the LE has been used as a probe 
to detect QPTs experimentally; at the same time, using 
a perturbative study in the short-time limit, the scaling 
relation T ~ (A)~ 22l/ valid close to a QCP (at A = 0) 
has been proposed. Here, v and z are associated corre- 
lation length and dynamical exponents, respectively [3[. 
In contrast to these studies where the coupling between 
the qubit and the environment is chosen to be weak, it 
has been shown that in the limit of strong coupling the 
envelope of the echo becomes independent of the cou- 
pling strength which may arise due to quantum phase 



transition in the surrounding [33J, [3J] . Moreover the LE 
and the decoherence of the central spin has been stud- 
ied when the environmental transverse Ising spin chain is 
quenched across the QCP by varying the transverse field 
linearly in time |35( . 

The central spin model we consider here, consists of 
a two level central spin S coupled to an environment E 
which is chosen to be a spin-1/2 XY spin chain with 
anisotropic interactions and subjected to a transverse 
field, described by the Hamiltonian 



H f 



N 



i=\ 



Vx« +1 + J y afaf +1 + hat], 



(1) 
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where erf (a = x,y,z) are Pauli spin matrices, h is the 
transverse field and N is the total number of spins in 
E. The spin chain (fl| is exactly solvable using Jordan- 
Wigncr mapping from spins to fermions [36h39| . The 
phase diagram is shown in the Fig. (fTJ); the transition 
from the ferromagnetically ordered phase to the param- 
agnetic phase driven by the transverse field h is called 
the Ising transition and the transition between two fer- 
romagnetically ordered phase, with magnetic ordering in 
the x direction (FM X ) and the y direction (FM y ), respec- 
tively, driven by the anisotropy parameter 7 = J x — J y , 
is called the anisotropic transition. The anisotropic tran- 
sition lines extend from h = — (J x + J y ) to h = J x + J y 
along the 7 = axis. Both Ising and anisotropic crit- 
ical lines meet at the multicritical points as shown in 
the figure. We shall exploit the exact solvability of the 
XY spin chain to calculate the LE close to these critical 
points. We also note at the outset that the spin chain 
is to be studied under periodic boundary condition and 
the wave vector k takes discrete values k = 2irm/N with 
m = 1,2,. ..N/2 and the lattice spacing is set equal to 
unity. 

Earlier studies focussed on the case when the coupling 
of the central spin to the environment is through the 
transverse field h |28l - l30l ] and explored the behavior of 
LE close to the Ising critical point. Motivation behind 
the present work is to explore the short-time behavior, 
collapse and revival of the LE around the anisotropic 
critical point (ACP) and especially the multicritical point 
(MCP) MCi of the phase diagram. To achieve this we 
evaluate the LE by coupling the qubit to the anisotropy 
term and also one of the interactions of the spin chain. 
Finally, we conjecture a generic scaling form that should 
be valid close to a QCP at least in the short-time limit. 

In the next section (Sec. II), we consider the case when 
qubit is coupled to the anisotropy term and the behavior 
of the LE is explored; in Sec. Ill, the qubit is coupled to 
one of the interaction terms (J x ). Finally in the conclud- 
ing section, we discuss our results and conjecture some 
generic scaling relations. 



FIG. 1: The phase diagram of the anisotropic XY model 
in a transverse field with Hamiltonian given by fl]) in the 
h/(J x + J y ) - j/(J x + Jy) plane, where 7 = J x — J y . The 
vertical bold lines denote Ising transitions from the ferromag- 
netic phase to the paramagnetic phase (PM), whereas the 
horizontal bold line stands for the anisotropic phase transi- 
tion between two ferromagnetic phases FM X and FM y . The 
multicritical points at J x = J y and h = ±1 are denoted by 
MCi and MC2, respectively. 



II. QUBIT COUPLED TO THE ANISOTROPY 
TERM OF ENVIRONMENT HAMILTONIAN 

In this section, it would be useful to rewrite the Hamil- 
tonian as 

1 N 
He = - g EK 1 + i)<°i + x + M " 7K°i!u + 2^1,(2) 

1=1 

with the choice J x + J y = 1, and the anisotropy parameter 
7 = J x — J y . Denoting the ground and excited states of 
the central spin by \g) and |e), respectively, the coupling 
of the system S to the environment E can be chosen as 



iY 



H. 



SE 



\e)(e\Y,[<r>i + i-vW + i], 



(3) 



where one assumes that the excited state of the qubit 
couples to all the spins of the environmental spin chain. 
The Hamiltonian of the composite system (S + E) is then 
given by 



H, 



H K +H. 



SE 



(4) 



The form of the interaction Hamiltonian chosen in 
Eq. ©, enables one to analytically calculate the behav- 
ior of the LE close to the anisotropic transition line and 
also the MCP. 

Let us assume the spin S to be initially in a pure 
state, \cj)(0))s = c g \g) +c e \e), (with coefficients satisfy- 

|2 , |„ I 



ln S \c g 



1) and the environment E be in the 



ground state denoted by |<p(0,7)).e; the total wave func- 
tion of the composite system at time t = can then be 
written in the direct product form 



so that the reduced Hamiltonian ((9j) gets modified to 



|*(0)) = |^(0)) 5 <8)|v(0,7)>i 



(5) 



One finds that the evolution of the XY spin chain splits 
into two branches (i) \(p(t,"f)) = exp(— iH( / y)t)\(p(0,'y)) 
and (ii)|^(t,7 + 5)) = exp(-i#(7 + 5)*)|<p(0,7)); this 
implies that \ip(t,-f)) evolves with the Hamiltonian @ 
with the anisotropy parameter 7 and \<p(t, 7 + <5)) evolves 
with the same Hamiltonian but the anisotropy parameter 
modified to 7 + 6. The total wave function at an instant 
t is then given by 

|*(t)) = c g \g)®\ V {t, 1 ))+c e \e)®\^{t, 1 + 5)). (6) 

One therefore finds the decay of the LE given by |28[ 

L (7,i)HMi, 7 Mt,7 + ( S))| 2 

= |(^(0,7)|exp(iiJ(7)i)cxp(-z J ff( 7 + J)i)|^(0,7)>| 2 

= M0,7)|expHH( 7 + <5)i)|^(0,7))| 2 , 

(7) 

where we have used the fact that the |<^(0)) is an eigen- 
state of the Hamiltonian ([2]) with anisotropy parameter 

7- 

The Hamiltonian © can be exactly solved by Jordan- 
Wigncr (JW) transformation followed by Bogoliubov 
transformations [3614391 ] and can be written in the form 
H{l + 5) = Ektkh + S)(A{A k - 1/2) and F( 7 ) = 
J2k £ k(l)(BlBk — 1/2) where A^'s and B k s are Bogoli- 
ubov fermionic operators and 



£fc(7 + 5) = y/(h + cos fc) 2 + {(7 + 6) sin fc} 2 ; 



(8) 



clearly £^(7) = £fc(7 + S) with 5 = 0. Here we have 
considered periodic boundary condition, the wave vector 
fc takes discrete values fc = 2-Km/N with m = 1,2, ...N/2 
( N is assumed to be even and also lattice spacing is set 
equal to one). 

In fact, under the JW transformation the Hamiltonian 
([2]) gets reduced to direct product of decoupled 2x2 
Hamiltonians for each momentum k which in the basis 
|0) (vacuum state) and |fc,— fc) (two JW fermion state) 
can written as 



H k {l) 



h + cos fc 17 sin fc 
— i7sinfc — (/i + cosfc) 



(9) 



In the current problem in which the LE is calculated as a 
function of 7, one makes resort to a basis transformation 
to |0) and |1), such that 



1 



V2 



(|0>+i|fc,-fc» 



(|0)-i|fc,-fc)), 



H k {n) 



7 sin fc h + cos fc 
h + cos fc —7 sin fc 



(10) 



The ground state of H(-f) and H (7 + 5) can be written 
in the form 

IW.O)) =cos— ^— |0)-zsm— — 11) 



lv(T + o, 0)) =cos 10) -ism 11), 

where tan0fc(7 + (5) = (/i + cosfc)/{(7 + 5) sinfc} and 
tan6> fc (7) = tan6» fc (7 + 6)\s=o- 

The Bogolioubov operators are related to the JW op- 
erators through the relation [36[ 



6 k ( 1 + S) . . 9^ + 5) t 
A}. = cos afe — 1 sin a 



-fc- 



(11) 



where a^-s are the Fourier transform of the JW opera- 
tors as derived through the JW transformations of spins. 
Using Eq. (jlip . one can further arrive at a relation con- 
necting the Bogolioubov operators 



B k = cos(a k )A k - i$\n(ak)A^_ k 



(12) 



where, a* = [0 fc ( 7 ) - fc ( 7 + <5)]/2 

Noting the fact that A k \tp(~/ + 6,0)) = and 
-Bfc|</?( 7 ,0)) = for all fc, one can use the Eq. (fT2")) 
to establish a connection between the ground states 
1 99(7 + 6,0)) and (93(7,0)) given by 

Mt,o)) = n[ cos K)+* sin K)4^ t -fc]^(o ) 7+^))-(i3) 

fc>0 

Substituting Eq. (|T3|) to the Eq. (JT]), we find the expres- 
sion for the LE given by 

L{l,t) = ]]L k = JJ[l-sin 2 (2a fc )sin 2 (£ fc ( 7 + ( 5)t)](14) 



fc>0 



fc>0 



We shall use Eq. ([14]) to calculate the LE as a function 
of parameters 7 and h, especially close the quantum crit- 
ical points. As shown in Fig. @, the LE as a function of 
7 (with the transverse field h < 1) exhibits a sharp dip 
near the anisotropic critical line (7 = 0). In contrary, 
when h = 1, and 7 is changed, we once again observe a 
sharp dip near 7 = which in this case happens to be the 
MCP MC\ (7 = 0, h = 1) as shown in the phase diagram 
Fig. (pj; changing 7 with h = 1 implies that we are in 
fact probing the behavior of LE along the gapless Ising 
critical line J4Q|. It should also be emphasized here that 
although the spin chain lies entirely on the critical line 
when the MCP is approached, we observe a substantial 
dip only at the MCP which suggests that the MCP is 
apparently playing the role of a dominant critical point 
in determining the temporal behavior of the LE [4l| . 
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FIG. 2: The LE as a function of 7 for h — 1; one observes a 
sharp dip around the MCP (7 = —5). Inset shows a similar 
dip around the aniosotropic critical line (with h = 0.5). 



the Hamiltonian if ((7 + S) = 0) is critical. We observe 
the collapse and revival of LE with time which is an 
indicator of quantum criticality as shown in fig. (|3|). 
It should be emphasized that when the size of the spin 
chain (E) is doubled keeping 6 fixed, the time period of 
collapse and revival also gets doubled; this confirms the 
scaling behavior mentioned above which is also observed 
at the Ising critical point |28j . 

The quasi period of oscillations can also be cal- 
culated in the following way. From Eq. (|14p . we 
find that the mode k = k c + 2-k/N gives dominant 
contribution for t — >• 00 so that for large N limit 
one can expand e\ in the form e\ = h + cosk ps 
\/I — h 2 2ir/N. We have also chosen 7 = —S such 
that $k(l + 6) = 7r/2 which makes 



sin 2a.k = 



7 2 sin k 



7 2 sin k + (h + cos k)' 



1. 



(17) 



Therefore from Eqs. (|14|) and (|17|). it is clear that 
oscillations in £(7, t) arises due to sin 2 e\t term pro- 
viding the time period 



A. Anisotropic Critical Point (ACP) 

The correlation length exponent v and the dynamical 
exponent z associated with the ACP is the same as those 
of the Ising transition, i.e., v = z = 1 and one therefore 
expects that the behavior of LE should be similar to that 
close to the Ising transition J28|. However, one needs to 
consider the fact that at the ACP the energy gap vanishes 
at 7 = for a critical mode k c = cos~ 1 (— h). 

As mentioned the decay of LE at short time is charac- 
terized by the critical exponents of the associated QCP. 
To calculate the short time behavior close to an ACP, we 
define a cutoff K c such that only modes up to this cutoff 
are incorporated in calculating the LE [28[ which is then 
given by L c (j, t) = Y\ k > £fe, and one defines 



S(j,t)=lnL c 



K c 

El 

fc>0 



In Li 



(15) 



we find sin 2 e^t : 



Expanding around the critical mode k c 
{ 1 + 8) 2 kH 2 and sin 2 (2a fc ) w fc 2 (5 2 /{ 7 2 ( 7 + 5) 2 } where 
we have relabeled k — k c as k; these lead to S("f,t) ps 
— J2k>o (k5At) 2 /j 2 . We therefore arrive at an exponen- 
tial decay of LE in the short time limit given by 



L c ( 7 ,i)«exp(-rt 2 ) 



(16) 



where, T = 5 2 E(K c )k 2 />y 2 and, E(K C ) 
{47r 2 A^ c (A^ c + l)(2A^ c + I)}/6Ar 2 (where N c is inte- 
ger nearest to NK c /2tt). From above equation (fT7j)) it 
is clear that in this case L c remains invariant under the 
transformation N — > Na, 8 — > S/a and t —$■ ta, with 
a being some integer. We now proceed to study the 
time evolution of LE with h = 0.5 and 7 = — 5 so that 



T = 



N 



2VT^h 2 



(18) 



This again shows that the time period of oscilla- 
tion of the LE is proportional to the size N of the 
environmental spin chain as shown in the Fig. ([3|). 
Eq. (|18p also shows that the time period diverges 
as h — >• I. This originates for the fact that for 
h = I, the spin chain lies on the gapless Ising crit- 
ical line, a situation which we are going to discuss 
in next sub-section. 

We note that the decay rate T scales as 7 -2 which is 
consistent with the scaling given in [32[ since zu = 1 for 
the transition across the anisotropic transition line also 

Si. 




FIG. 3: The variation of the LE as a function of time at 
the ACP for 7 = -5, h = 0.5 and 5 = 0.01. The collapse 
and revival of LE is indicator of a QPT and the quasiperiod 
of the LE is proportional to the size N of the environmental 
spin chain 



B. Multicritical Point (MCP) 

As mentioned already, we set h = 1, and ap- 
proach the MCP by changing 7 along the Ising crit- 
ical line. Expanding sme^t and sin(2afc) near MCP 
around k = n; sin 2 e k e t sa (7 + 5) 2 k 2 t 2 and, sin 2 (2a k ) ~ 
fc 2 <5 2 /{47 2 (7 + S) 2 }; k — > (fc— 7r), one finds the short-time 
decay of LE given by 



L c ( 7 ,t)«exp(-rt 2 ) 



(19) 



where, T = S 2 E{K c )/4:j 2 and, £(if c ) 
{(1/5)A C 5 + (1/2)A C 4 + (1/3)A C 3 - (1/30)A C }/7V 4 
(where N c is integer nearest to NK c /2ir). Equa- 
tion p^|) helps in providing analytical scaling for LE 
i.e., L c is invariant under transformation N — > Na, 
S —> 5 /a 2 and t — > to? . This scaling has to be contrasted 
with the scaling of L c close to the ACP presented in 
the previous section. We note that at the MCP, the 
minimum energy gap scales as (fc — n) 2 so that z = 2 
whereas near an ACP , it scales linearly as (k — k c ) 
with 2 = 1. This difference in the dynamical exponent 
is the reason behind different scaling observed in the 
short-time limit. The collapse and revival of LE as a 
function of time is shown in Fig. Q for different system 
sizes and fixed 5; this confirms the scaling observed in 
the short-time limit. At the same time, we note that 
r ~ 7~ 2 as the exponent zv = 1, even for transition 
across the MCP. 

To calculate the time-period of oscillation, we 
again proceed using the same line of arguments 
given in section A. At the MCP (h=l), e k e = 
h + cosfc Ks 2ir 2 /N 2 and similarly for 7 = —S, it can 
be shown that sin 2 2a k ~ 1. The time period of os- 
cillations in L(j, t) is therefore given by T ss N 2 /2w, 
which confirms that the LE oscillates with period 
is proportional to N 2 at MCP (see Fig. (g])). 
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FIG. 4: The LE is shown as a function of time at the MCP 
j — —5 , h — 1 and 8 = 0.01. The quasiperiod of the collapse 
and revival is proportional to iV 2 according to the scaling 
relation N — s> iVa, 8 — > 8/ a 2 and t — ► ta 2 , discussed in the 
text. We note that the collapse and revival close to the MCP 
is not a smooth function of time. 



III. QUBIT COUPLED TO THE INTERACTION 

TERM OF THE ENVIRONMENT 

HAMILTONIAN 

In this section we shall choose the form of the XY 
Hamiltonian given in Eq. ( fT|; t ransforming to a new state 
of basis vectors defined by [42| 



\e lk ) =sin(k/2)\0) + icos(k/2)\k,-k) 



e 2k 



cos(fc/2)|0) -isin(k/2)\k,-k), 



one can rewrite the reduce 2x2 Hamiltonian © Hk in 
the form 



Jx + Jy cos 2k + h cos k 
J y sin2fc + /isinfc 



J y sin2fc + /isinfc 



-(J x + Jy cos 2k + h cos fc) 
We choose the coupling term given by 



.(20) 



A r 



Hse = S\e){e\Y^[«+i], 



(21) 



;=i 



Therefore total Hamiltonian (spin and environment) be- 
comes 

N 

H e = - £[(•/, + S\e)(e\)afaf +1 + J y o\a\ +1 + haf], (22) 



The advantage of selecting such a coupling is that it 
enables us to explore the MCP and ACP via different 
paths and compare the results with the previous case, 
e.g., if one chooses J x = 2h y , the MCP is approached 
along a linear path when J x is changed unlike the pre- 
vious case when it is approached along the Ising critical 
line. 

Following identical mathematical steps as described in 
the previous section, one can find that the expression of 
the LE is given by 

L(Jx,t) = l[L k = [][l-sin 2 (2a fe )sin 2 (e fc (J ;c + ( 5)t)](23) 
k>0 k>a 

where, a k = [6 k (J x ) - 6 k (J x + S)}/2, and 

J y sin(2fc) + h sin(fc) 



1k(Jx+S) = arctan[ 



J x + 5 + J y cos(2fc) + /icos(fc) 



(24) 



The energy spectrum given in Eq. ((5J can be rewritten 
as 

£k(Jx + 5) = [(Jy sin 2fc + /isinfc) 2 

+ (Jx +5 + J y cos2k + h cos k) 2 } 1/2 . (25) 

Let us first explore the LE close to different critical 
points; refereeing to fig §5$ and ©, we find that there 
is a sharp dip in LE wherever the parameters values are 
such that the system is close to a critical point. For 
example, in Fig. ([5]), we have varied J x keeping h and J y 
fixed and h < 2J y (= 2) such that we observe dips at two 
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FIG. 5: The LE plotted as a function of J x shows dips 
around the Ising critical points (J x = —0.2, —1.8) as well as 
the ACP (J x = 1) with J v = 1 and h = 0.8(< 2J y ). Inset 
shows that for h = 2.2(> 2J y ), there are dips only at the Ising 
critical points as the variation of J x does not take the spin 
chain across the anisotropic critical line. 
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FIG. 6: The interaction J x is varied with h fixed to h = 2J y = 
2. The LE shows a dip at Ising critical point(J :E = —3) and 
also at the MCP (J x = 1). 



Ising critical points and also at the anisotropic critical 
point; for h > 2J y , in contrast, one observes dips only 
at the Ising critical points as the anisotropic transition 
point is not crossed in the process of changing J x . For 
h = 2J yi one observes dips at the Ising critical point and 
the MCP as J x is varied (Fig. (|6])). Equipped with these 
observations, we now proceed to study the short time 
decay of LE close to these critical points. 



shown in Figs. (0 and <j8j> . This confirms that the scaling 
of LE does not depend on how the central spin is cou- 
pled to the environment rather it is hypersensitive to the 
proximity to a critical point of the environment. 
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FIG. 7: The collapse and revival of the LE at the ACP 



(./, 



d , h = 0.8, J v 



1). The inset shows the same 



behaviour at the Ising critical point (h = 2.2, J x 
J v — 1). In both the cases 8 = 0.01 
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FIG. 8: The time variation of LE at MCP ( J x = 1-5, h = 2 = 
2J y ) is shown . The quasiperiod of LE is again proportional 
to iV 2 as reported in Sec. II. 



B. Close to the MCP 



A. Short time behavior 

Near the critical points and the MCP, we have identi- 
cal short-time behavior and scaling with respect to N, S 
and t as already reported in the previous section. These 
are corroborated by numerical estimation of collapse and 
revival close to the critical and the multicritical point as 



It is well known that for a finite XY spin chain, 
there exist quasicritical points on the ferromagnetic side 
close to the MCP; the energy gap is locally minimum 
at these quasicritical points and it scales k 3 in contrast 
to the scaling k 2 at the MCP. In the limit of N — > oo, 
all these quasicritical points approach the MCP. These 
quasicritical points and exponents associated with them 
have been found to dictate the scaling of the defect den- 



sity following a slow quench across the MCP [4l|, |43| 
and also the scaling of fidelity susceptibility close to 
it (4J|. We shall now explore the collapse and revival 
of the LE fixing the parameters such that J x + 5 is 
right at a quasicritical point. For modes k ra n, one 
can use the simplification, sin 2 e\t w (J x + S — J y ) 2 t 2 
and, sin 2 (2a fc ) w 4J 2 fc 6 <5 2 /{(J ;c - J^V* + <J - J y ) 2 }. 
We therefore get a similar exponential decay of the LE 
L c {J x ,t) w cxp(-R 2 ) with 



r = 



4J 2 <5 2 



s(i<r c 



(J. - Jy) 2 



and E(K C ) = 



A(N C ) 



(26) 



where A(N C ) = {l/7)Nj+(l/2)N^+{l/2)N^~{l/6)N 3 + 
(l/42)A r c , and as defined previously, iV c is integer nearest 
to NK c /2tt. The above equation (f2"6")l , shows a very in- 
teresting scaling behavior of the LE TV — > Net, 5 —> S/a 3 
and t — > to 3 which is different from the scaling ob- 
served at the MCP. Simillar to previous cases, at 
the quasicritical point e\ = h + cosk ps 16ir 3 /3N 3 
with h = 2 J v = 2; for J x = 1 - S + 4tt 2 /N 2 and large 
N, it can be easily shown that sin 2 2au ~ 1. The 
time period of oscillations in L(J x ,t) is therefore 
given by T « iV 3 /167r 2 , which verifies the fact that 
LE oscillates with period proportional to A^ 3 at 
quasi critical point (as shown in Fig. (|9j) ) . 

The collapse and revival of the LE as a function of 
time supports the scaling behavior analytically obtained 
in the short time limit (see Fig. ([9])). Comparing Eq. (|26|) 
with the form of decay rate T given in Eq. (jf6|) . we find 
that in both the cases T ~ I/7 2 ; this is because at a 
quasicritical point one can define an effective dynamical 



exponent z„ 



3 v a 



1/3 such that v q 



1 41, 4J 



Moreover, we find that the quasipcriod scales as N z 
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FIG. 9: Collapse and revival of the LE at a quasicritical 
point (J x = 1 — 5 + Ait 2 /N 2 , h = 2 and J y — 1) as denned in 
the text. The quasiperiod scales with iV 3 in contrast to ~ iV 2 



at the MCP. This is consistent with the scaling N 

5 -¥ S/a 3 and i -^ to 3 . 



Na, 



IV. CONCLUSION 

In this paper, a spin-1/2 (qubit) is coupled to the en- 
vironment which is chosen to be a spin-1/2 XY spin 
chain and the temporal behavior of the LE is studied. 
The coupling is done in such a way that enables us to 
study the LE close to the ACP as well as the MCP of 
the phase diagram and these points are approached in 
different fashions, e.g., the MCP is approached along the 
Ising critical line in Sec. II while in Sec. Ill it is ap- 
proached following a linear path. We find that close to 
the ACP, the evolution of the LE is identical to that re- 
ported in the ref. [28|]. However, around the MCP, we 
observe that the quasiperiod of the collapse and revival 
of the LE as a function of time scales as N 2 where N is 
the size of the environmental spin chain. We attribute 
this to the fact that the dynamical exponent z associ- 
ated with the MCP is two. To justify this conjecture, 
we have estimated the scaling of the decay rate T and 
also the period of the collapse and revival of the LE at 
a quasi-critical point on the ferromagnetic side of the 
MCP. We find that quasiperiod scales as A^ 3 . It should 
be noted here that at the quasicritical point, the mini- 
mum gap scales with the system size as A^ 3 and hence 
one can define an equivalent dynamical exponent z qc = 3 
[43j . In Sec. II, even though the MCP is approached 
along a gapless critical line, a sharp dip in the LE is ob- 
served only around the MCP where the decay of energy 
gap with the system size is faster (~ 1/N 2 ) with respect 
to that near the Ising or anisotropic critical point. We 
observe that the collapse and revival of the LE at MCP is 
not a smooth function of time which is attributed to the 
fact that in Sec II the spin chain is always close to the 
Ising critical line whereas in Sec. Ill quasicritical points 
are likely to influence the temporal evolution of the LE. 
These quasi-critical points, on the other hand, are ex- 
pected to be related to the proximity to the critical line 
of the finite-momentum anisotropic transition. 

Although wc have studied an intcgrablc spin chain re- 
ducible to direct product of two-level system, our studies 
indicate the possibility of some interesting scaling behav- 
ior. We see that in all the cases studied here, the LE de- 
cays exponentially close to the critical point in the short 
time limit with the decay rate T scaling as T ~ \~ 2zv i.e., 
our studies support the scaling proposed in [32| based on 
perturbative calculations and a Landau-Zcner argument. 
Moreover, we find the quasiperiod of the collapse and 
revival of the LE at the critical point scales as N z ; we 
note that the dynamical exponent z determines how does 
the minimum energy gap vanishes with increasing system 
size (~ A^ -2 ) at the QCP. At a quasicritical point the ef- 
fective dynamical exponent z qc = 3 is found to determine 
the scaling of the quasiperiod of collapse and revival with 
the system size. 

Finally, we comment on the decoherence of the central 
spin during time evolution which is calculated using its 
reduced density matrix [28[. The off-diagonal terms of 
the reduced density matrix is given by c*c e d(t) and its 



hermitian conjugate where the decoherence factor d(t) is 
connected to the LE through the relation L(t) = \d(t)\ 2 
[35j | . The vanishing of the LE around to the QCP there- 
fore implies a complete loss of coherence and therefore 
the qubit makes transition to a mixed state even though 
initial state is chosen to be pure. On the other hand, 
away from the QCP LE stays close to unity, thus the 
purity of the qubit state is retained. Our studies reveal 
that close to the MCP, T ~ 1/N 2 in the short time limit, 
implying a faster loss of coherence with the increasing 
system size when the environment E is close to a MCP 



than when it is close to a QCP. The loss is even faster 
when the spin chain sits at a quasicritical point close to 
the MCP. This faster loss of coherence with the system 
size, we believe, is a note-worthy observation. 
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